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ABSTRACT 

We present a 3D general circulation model of Pluto and Triton's at- 
mospheres, which uses radiative-conductive-convective forcing. In both the 
Pluto and Triton models, an easterly (prograde) jet is present at the equator 
with a maximum magnitude of 10-12 m s _1 and 4 m s _1 , respectively. Neither 
atmosphere shows any significant overturning circulation in the meridional 
and vertical directions. Rather, it is horizontal motions (mean circulation 
and transient waves) that transport heat meridionally at a magnitude of 
1 and 3 xlO 7 W at Pluto's autumn equinox and winter solstice, respec- 
tively (seasons referenced to the Northern Hemisphere). The meridional and 
dayside-nightside temperature contrast is small (<5 K). We find that the 
lack of vertical motion can be explained on Pluto by the strong temperature 
inversion in the lower atmosphere. The height of the Voyager 2 plumes on 
Triton can be explained by the dynamical properties of the lower atmosphere 
alone (i.e., strong wind shear) and does not require a thermally defined tro- 
posphere (i.e., temperature decreasing with height at the surface underlying 
a region of temperature increasing with height). The model results are com- 
pared with Pluto stellar occultation light curve data from 1988, 2002, 2006, 
and 2007 and Triton light curve data from 1997. 

Keywords: ATMOSPHERES, DYNAMICS ; OCCULTATIONS ; PLUTO, 
ATMOSPHERE ; ATMOSPHERES, STRUCTURE ; TRITON 
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1. Introduction 



Pluto and Triton are sometimes referred to as "sister" worlds. They have 
similar radii of 1132-1200 km for Pluto (Young and Binzel, 1994; Tholen and 
Buie, 1990, 1997; Buratti et al., 1995; Reinsch and Fes tou, 1994; Zalucha et 



al. 2011a,b) and 1353 km for Triton ( iThomasl . 120001 ). Their atmospheric 



compositions are primarily N2 with trace amounts of CH4 and CO (and 



CQ 9 for Triton) (IQwen et all Il993l: lYoung et al 



Herbert and Sandel 



1991 



Cruikshank et al. 



1993 



997 



Greaves et al. 



Lellouch et al 



2009 



2011 



2010 



201 ll ). Both have a rotation rate of about 6 Earth days and mean surface 
pressure of tens of microbars (Broadfoot et al., 1989; Lellouch et al., 2009; 
Zalucha and Gulbis, 2012). Additionally, both receive weak insolation due 
to their large distance from the Sun. Thus, their atmospheres are in very 
similar regimes (unlike Mars and Venus, or Earth and Titan), and we may 
apply the same model framework to both of these bodies. 

A large amount of analysis and modeling has been done on the ID ver- 
tical thermal structure of Pluto and Triton's atmospheres and surfaces (e.g., 
Hansen and Paige, 1992,1996; Yelle and Lunine, 1989; Hubbard et al., 1990; 
Stansberry et al., 1994; Strobel et al., 1996; Elliot and Young, 1992; Elliot et 
al., 2003; Zalucha et al., 2011a,b). While it is mathematically possible that 
these atmospheres' vertical temperature structure is set by energetic pro- 
cesses alone, it is rather unlikely that Pluto and Triton have no significant 
atmospheric circulation based on observations and physical theories of plan- 
etary atmospheres. Though somewhat inhibited by the relative difficulty in 



observing atmospheric circulation, little work has been done exploring the at- 
mospheric circulation of Pluto and Triton. Occultation studies of Pluto have 



provid ed some measure of t he wind speed (IPerson et all 



waves ( IHubbard et al. 



20081 ) and gravity 



20091 ) on Pluto, while the Voyager 2 flyb y of Triton re 



suite d in observations of surface-based vertical jets of aero sols (ISmith et al 



1989) that were sheared off aloft by high altitude winds (ISoderblom et al. 



19901 ) . Subsequent ana 



surfac e wind streaks by 



ysis was performed on the Voyager 2 obser vations of 



Hansen et al. 



(1990) 



(jl990l ) and plume dynamics by 



Ingersoll 



Recently, numerically sophisticated general circulation models (GCMs) 
have been applie d to Pluto (Zalucha and Gulbis, 2012; Michaels and Young, 



Miller et al. 



2011) and Tr iton ( Mueller- Wodarg et al. 



2001 



Vangvichith and Forget 



2011 



20111 ) . GCMs not only include the physics of fluid motion, but 
are also coupled with surface processes and thermodynamical processes such 
as radiation, conduction, and condensation. Despite the sudden proliferation 
of Pluto GCMs (PGCMs) and Triton GCMs (TrGCMs), all currently exist- 
ing ones use different physical assumptions and a consensus about Pluto's 
atmospheric circulation is far from being reached. Some of the difficulty 



may be attributed to Pluto and Triton's steep, 



version ( Yelle and Lunin 



1989 



Hubbard et al. 



qw-alt i tude temperature in 



19901 : 



Strobel et al 



19961 ). 



which makes analogies between other Solar System bodies such as Mars and 
Titan difficult beyond the very simplest arguments. 

This paper continues the work of Zalucha and Gulbis (2012, hereinafter 
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referred to as ZG12) from a 2D PGCM model to a 3D PGCM and also applies 
this particular model to Triton for the first time. The ZG12 PGCM/TrGCM 
uses t he dynamical core of t he Massachusetts Institute of Technology ( MIT) 



GCM (Marshall et al 



19971 ) and radiative-conductive forcing scheme of lYelle and Lunine 



(119891 ) (now augmented with the convective forcing of Zalucha et al., 2011b). 
The PGCM/TrGCM makes implicit assumptions about the surface, subsur- 
face, and condensation cycle (elaborated in ZG12). Surface and subsurface 
parameters (such as surface emissivity, surface albedo, surface thermal iner- 
tia, and surface frost history) are not well constrained and represent a large 
area of parameter space that is not possible to analyze in a single paper and 
is left for future work. 

Stellar occultation observations represent a significant fraction of the 
available data for Pluto, and for Triton to a lesser extent. We use the for- 
ward method for validating the PGCM and TrGCM with the data; namely, 
we take pressure-temperature profiles from the PGCM/TrGCM, calculate a 
model light curve, and assess any differences in normalized flux vs. radius (or 



time) space. Temperature profiles derived from inverse methods (lEUiot et al 



20031 ) may contain systematic errors due to the assumption of an 



upper 



19921 ) must 



boundary condition, while idealized methods (lElliot and Young , 
make assumptions about the variation of temperature with height. The 
GCM to light curve conversion has been used numerous times (Zalucha et 
al., 2011a,b; ZG12) and we are increasingly confident of this technique. 
In Section [2] we describe our PGCM and TrGCM setup and briefly sum- 
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marize the light curve model. In Section EJ we compare the PGCM and 
TrGCM output with stellar occultation light curves and present the corre- 
sponding temperature and zonal wind fields. We also show the meridional 
heat transport on Pluto. In Section @] we discuss the effect of Pluto's lapse 
rate on atmospheric circulation and provide evidence that Triton's lower at- 
mospheric structure is dynamically (rather than thermally) maintained. 

2. Model setup 

2.1. General circulation model 

Both the PGCM and the TrGCM use the same dynamical core, based 
on the MIT GCM. The configuration of the dynamical core was described 
in ZG12 and will only be briefly summarized here. The MIT GCM solves 
the primitive equations of geophysical fluid dynamics (i.e., the Navier-Stokes 
equations) under the assumptions of conservation of mass, the ideal gas law, 
and conservation of energy in the presence of an external source and sink 
term) on a sphere using the finite volume method. The model is hydro- 
static in the vertical and compressible. The vertical grid is based on an 7] 
coordinate and the surface topography is flat. Boundary layer friction is 
represented by a simple drag law (linearly dependent on the horizontal ve- 
locity) that decreases with height, reaches zero at the top of the boundary 
layer, and is zero at all levels above. The composition of the atmosphere is 
assumed to be primarily N2, with trace amounts of CH4 that is radiatively 
active but not advected. In this version of the model, mass is not permitted 
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to be exchanged between the surface and the atmosphere (but will be in fu- 
ture versions). When the atmospheric temperature is diagnosed to fall below 
the N 2 freezing temperature, it is instantaneously reset back to the freezing 
temperature. Surface thermal inertia and albedo are globally and tempo- 
rally constant (this assumption will also be relaxed in the future). Surface 
temperature is held at the N 2 freezing tempe rature, consistent with observa 



tions that Pluto and Triton arc frost-covered (Owen ct al. 



1993 



Tryka et al. 



19931 ). The implications of this assumption are discussed in Section I3TT1 The 
latitude and longitude coordinates we use to describe locations on Pluto and 
Triton's surfaces follow the convention that Pluto's north pole is in the same 
hemisphere as the ecliptic north pole with longitude is increasing eastward. 
This convention is equivalent to stating that both Pluto and Triton rotate in 
a retrograde fashion (i.e., clockwise when viewed from north). 

An important improvement since ZG12 is the expansion of the model 
domain to three spat ial dimensions. In the horiz ontal, the MIT GCM uses 



a cubed-sphere grid (lAdcroft and Campin 



20041 ) with 32 x 32 points per 



cube face, equivalent to a grid spacing of 2.8° at the equator. Compared to 
the more common cylindrical projection grid (i.e., a latitude/longitude grid), 
this type of horizontal grid eliminates singularities at the poles that force 
meridional winds to zero and removes the requirement for Fourier filtering 
in the high latitudes (in order to maintain a practical time step). 



As in ZG12, we use the radiative-conductive scheme of 



Yelle and Lunine 



( 119891 ) to calculate the external heating and cooling terms due to CH 4 and 
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molecular conduction (now including the longitudinal dependence in inso- 
lation). This scheme captures non-local thermodynamic equilibrium (non- 
LTE) effects and the sharp stratospheric temperature inversion, while being 



computationally practical in a GCM. 



existence of a troposphere on Tr iton (ISmith et al 



Voyager 2 o 



Dservations indicate the 



ified the 



19891 ). so we have mod- 



Yelle and Lunind ( 119891 ) model to include this effect following the 



procedure described in Zalucha et al. (2011b). The heat balance equation is 



dT d ( dT\ 



(i) 



where p is atmosperic density, c p is the specific heat at constant pressure, T is 
temperature, z is altitude, the third term on the RHS is the conduction term 
with thermal conductivity coefficient k a and thermal conductivity exponent 
a, Q is the non-LTE heatin g by CH 4 at 3.3 /xm, and I is the non-LTE cool- 



ing by CH4 at 7.8 /im. The 



Yelle and Lunind (119891 ) model used to specify 



Q and I doe s not include CO, which has been recent 



atmosphere (IGreaves et al. 



2011 



Lellouch et al. 



y observed in Pluto's 



20111). It is difficult to an- 



alytically express the cooling effects of CO; thus CO has not been included 



at this stage owi ng to computa t iona 



of CO using the 
work. 



Strobel et al. 



simplicity. We plan to add the effects 



(|l996l ) radiative-conductive model in future 



In Eq. [TJ Cdiff is the term that prescribes the troposphere via convection 
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or 

where 9 e is the equivalent potential temperature (due to condensation or 
sublimation of N 2 ) and Kc is the convection diffusion coefficient. Note that 
here we have used equivalent potential temperature assuming implicit con- 
densation and sublimation of N 2 in order to mathematically represent the 
troposphere. As stated previously, no explicit exchange of mass occurs in 
the continuity equation (in the model's dynamical core). Assuming that the 
tropospheric temperature is primarily determined by convection, the vertical 
temperature profile will be equal to the moist (with respect to N 2 ) adiabatic 
lapse rate of —dT/dz = 0.09 K km -1 . Specifically, 



K c = 20 



1 — tanh 



z — h, 



(3) 



where z has units of km and h c is the troposphere critical height (in km), 
interpreted as the altitude to which convective effects are significant. h c 
varies monotonically with the depth of the troposphere (i.e., larger h c corre- 
sponds to larger troposphere depth), but the two are not necessarily equal. 
The functional form of Eq. [3] was chosen based on the structure of Earth's 
troposphere as in Zalucha et al. (2011b). 

It has been known for many decades that terrestrial climate models need 
a momentum sink in the mesosphere to match observations . Initially, this 



was parameterized with Rayleigh friction (e.g., iLeovyl . 



19641 ). a momentum 
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sink that was linearly dependent on velocity and increased in magnitude with 
height. More elaborate schemes varied this forcing in time and space until the 
desired result was achieved. Eventually, it was discovered that the breaking 
of subgrid-scale gravity waves (also known as buoyancy or inte rnal waves' 



was the physical me c hanism 



1978 



Lindzenl . 



1981 



Holton , 



responsible 



1982 



1983 



br the momentu m drag (IHoughton . 



Matsuntol . 



19821 ). While the theory 



behind linear, small-amplitude, idealized gravity waves is well-understood, 
these waves present a problem for global climate models because their typical 
wavelength is smaller than the grid spacing. Hence, the amount of momen- 
tum deposition due to breaking waves must be parameterized. 

V arious gravity wave parameterizations hav e been imple mented for Earth 



Lindzen and Holton 



1968 



Lindzen 



198lh . Mars (e.g.. 



Medvedev et al. 



201 lh . and Venus (e.g., Zalucha et al., 2013) with varying degrees of success. 
In all cases, the parameterizations must be tuned to observations (usually 
of perturbations in vertical profiles of temperature or density). Observa- 
tions of such waves on Pluto and Triton are almost nonexistent (but see 



Hubbard et al 



20091 for one example). Thus, in the PGCM and TrGCM 



we must resort to Rayleigh friction. For numerical stability, a surprisingly 
large momentum sink is required in the uppermost five levels of the model 
domain. In order of decreasing height, the drag coefficients are 0.02083 s _1 , 
0.00694 s _1 , 0.00231 s~\ 0.00072 s" 1 , and 0.00024 s" 1 ; for comparison, typ- 
ical values for Mars GCMs in the uppermost three levels are 10.1 x 10 -5 s _1 , 
3.4 x 10~ 5 s -1 , and 1.1 x 10~ 5 s -1 . The 2D PGCM did not require such 
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strong friction (in fact it was stable even with no Rayleigh friction). The 
reason for such strong Rayleigh friction might indicate that vertically propa- 
gating gravity waves deposit relatively large amounts of momentum in Pluto 
and Triton's upper atmosphere. Were we able to incorporate a true grav- 
ity wave parametrization (i.e., if we had observational information about 
the wave structure on Pluton and Triton), the strength of the drag would 
potentially depend on longitude and latitude. 

In practice, the Rayleigh friction layer (often referred to as the sponge 
layer) is not treated as giving physically meaningful results, and the model 
output in the sponge layer is ignored. Likewise, we do not present results 
from this layer in our results (Section E]) . The effect of choosing the top five 
layers as the sponge layer, rather than the top four or six, simply determines 
the altitude to which we are able to ascertain the circulations of Pluto and 
Triton's atmosphere. In the 2D version of the PGCM, a thicker sponge layer 
led to weaker zonal winds. We have not performed a detailed analysis on 
the effect of the strength of the momentum sink in either the 2D or the 3D 
PGCM. As stated above, the 2D PGCM was stable without a sponge layer at 
all, which is due to the fact that zonally propagating waves are not possible 
in a 2D GCM. 

Both the PGCM and TrGCM were run for a 15 Earth- year spin up period 
(as explained in ZG12), then continued through the year 2007. Since the 
spin up time is determined primarily by the radiative-conductive relaxation 
rate, the addition of a third spatial dimension does not change the spin 
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up ti me. The first occu ltation dataset for Pluto was obtained in the year 



1988 (Elliot et all . Il989f ). so the PGCM simulations begin in the year 1973. 



Fewer occultation datasets exist for Tri ton; the hig 



r est sig nal-to-noise ratio 



(SNR) dataset was obtained in 1997 (jEUiot et al 



19981 ). so the TrGCM 



simulations begin in the year 1982. In both cases, the atmospheres were 
initialized in radiative-conductive-convective equilibrium (i.e., steady state 
solution to the radiative-conductive-convective scheme with no wind) and 
with globally constant surface pressure. 

Table [TJ shows the bulk parameters of the PGCM and TrGCM. For Pluto, 
we performed a parameter sweep of surface pressure in the range of 4-28 /ibar 
in intervals of 4 /zbar. The CH 4 mixing ratio was set at 0.2%, 0.6%, and 1%. 
The ranges and intervals we re chosen in accorden ce with previous analyses 



and observations (e.g., ZG12, 



Lellouch et al. 



20091 ). We determined the best- 



fit surface pressure/ CH4 mixing ratio parameter combinations for a particular 
observational dataset by calculating a model stellar occultation light curve 
from each PGCM simulation (see Section 12 .21) and then finding the model 
run that has the minimum reduced \ 2 compared to the data. The existence 



of a troposphere on Pluto is debated. 



Stansberry et al. 



(119941 ) showed that 



the addition of a troposphere caused a kink in the l ight c urve (similar to 



that observed in the 1988 occultation by 



Elliot et al. 



19891 ). The flux level 
of the kink and depth of the light curv e could be changed by v arying the 



depth and lapse rate of the troposphere. 



Stansberry et al. 



( 119941 ) show that 



a troposphere was consistent with the data, but not required. Moreover, 
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Zalucha et al. (2011b) find better agreement for a model with no troposphere. 
Thus, we assume no troposphere in the PGCM. 

Triton was visited by Voyager 2 in 1989, so we used measured values 
from this encounter, along with more recent measurements, to configure the 
TrGCM. We then compare the results to the 4 Novembe r 1997 stellar occu l- 



tation observed by the Hubble Space Telescope (HST) (jElliot et al.l . Il998l ). 
which is the only Triton occultation dataset with the minimum SNR thresh- 
old to be used in our analysis (see Zalucha et al., 2011b for a definition of 
the minimum SNR). 

During the Voyager 2 encounter w ith Triton in 198 9 the a tmospheric pres- 



sure was determined to be 14 /zbar ( jBroadfoot et al. 



1989 ). The CH 4 par 



tial pressure was det e rmine d to be 2.45 nbar ([Herbert and Sande 



1991 



Strobel and Summers 



1995|). 



Elliot et al. 



(119981 ) observed from a stellar oc- 



cultation that occurred on 4 November 1997 that the pressure at the half 
light radius had increased significantly since the Voyager 2 encounter and 
measurements from a 1995 stellar occultation. They concluded that the sur- 
face pressure of Triton doubles every 10 years du ring the 1989-1997 pe riod of 



Lellouch et al. 



(120101 ) found 



observation. Using spectroscopic measurements 
that in 2009 the partial pressure of CH4 had increased to 9.8 ±3.7 nbars. 
Combining all this information and performing a linear interpolation, we in- 
fer that the surface pressure in 1997 was 25.2 /ibar and the mixing ratio 
of CH4 was 0.02%. Voyager 2 also observed plumes rising fro m the surface 



and then being sheared off at an altitude of approximately 8 km 



Smith et al. 



14 



(119891 ). T his has been inter preted as the depth of a thermally-defined tro- 



posphere 



Yelle et al. 



19951 ). Since the tropospheric depth has not been 
observed since the Voyager 2 encounter, we take the troposphere depth in 
the TrGCM to be 8 km. We reiterate that our method is not simply a reanal- 



ysis of Voyager 2 derived results. 



measurements (lElliot et al. 



1998 



jut synthesizes it alo ng with more recent 



Lellouch et al. 



20101). 



2.2. Light curve model 

The PGCM and TrGCM output pressure-temperature vertical profiles 
as a function of latitude, longitude, and time. Using hydrostatic balance, 
the ideal gas law, and the linear dependence between number density and 
refractivity, the refractivity v as a function of the body's radius r may be 
ob tained. Then the bending angle 9 and its vertical derivative d9/dr is given 



by 



Chamberlain and Elliot 



(119971 ): 



9{r) 



r dv{r') 
r' dr' 



dx 



(4) 



and 



d9(r) 
dr 



x 2 dvir') r 2 d 2 v{r') 



) 3 dr' 



dx. 



(5) 



where x 2 = r' 2 — r 2 is the path along the ray and r' is a dummy variable of 
integration. The normalized light curve flux is given by (assuming no line of 
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sight extinction) 



^ \l + D[d9(r)/dr]\\l + D6(r)/r\' ^ 



where D is the distance between the observer and body. In this model, we 
have assumed no line of sight extinction along the incoming stellar light ray. 
See Zalucha et al. (2011a) for a further explanation of this method. 

The forward method of light cur ye modeling is employed, where the GCM 



results are used as input into the 



Chamberlain and Elliot 



(119971 ) scheme, 



yielding a model, GCM-based light curve. The agreement between data and 
model is then evaluated in flux vs. radius (or flux vs. time) space. Idealized 
light curve modeling techniques provide temperature as a function of radius 



from light curve data. However, these techniques (such as lElliot and Young . 



19921 ) must make an assumption about the temperature structure (in this 
case that is follows a power law), which is usually extremely simplified for 
mathematical convenience and thus not accurate for planetary atmospheres 



in ge neral. Other techniques such as the inversion method (lElliot et al. 



20031 ) require an upper boundary condition, which is usually obtained from an 
idealized model and can introduce large systematic errors. Forward modeling 
can utilize atmospheric fields having any vertical structure and is therefore 
potentially more accurate. 
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3. Results 



3.1. Best- fit PGCM light curves 

We validate our PGCM using the following datasets, which are selected 
based on the criterion that they have the highest SNR per scale height (see 



Zalucha et al., 201 



(KAO) (Elliot et al. 



a): 9 June 1988 on the Kuiper Airborne Observatory 



19891 ) 0.9 m telescope; 21 August 2 002 at Mauna Kea 



on the Universi t y of Hawaii 2.2 m telescope (UH 2.2m) ( lEUiot et al 



2003 



Pasachoff et al 



20051 ); 12 June 2006 at Siding Spr ing, Australia on t he 2.3 m 



Australian National University (ANU) tele scope ( 



the 4 m Anglo- Australian Telescope (AAT) (lYoung et al. 



Elliot et al. 



Qlkin et al. 



20071 ) and on 



2008): and 31 June 



20091 ). 



2007 at Mt. John Observatory in New Zealand (J< 

In our current analysis we also performed PGCM simulations in which 
the surface temperature was either free to evolve (i.e., determined by the 
insolation on an non- volatile surface) or fixed (i.e., held fixed at the N 2 freez- 
ing temperature). The parameters of best-fit (i.e., minimum reduced y 2 ) 
are shown in Tables |2] and [3j respectively. Note that these values assume 
a surface radius of 1152 km, discussed below. The formal error bars (see 
ZG12 for details about this calculation) are smaller than the interval of the 
parameter sweep, so we take the interval of the parameter sweep to be the 
reported error bars. The fixed surface temperature simulations always have a 
lower reduced x 2 ; indicating that they are better fits. Moreover, the derived 
surface pressures for these runs agree better with previous analyses (Lellouch 
et al., 2009; Zalucha et al, 2011b). 
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To put into perspective the differences in the reduced x 2 values of the fixed 
vs. free surface temperature simulations, we point out that in all years, the 
reduced x 2 value for adjacent simulations in parameter space (i.e., the one- 
sigma error bars) for the fixed surface temperature simulations was still less 
than the free surface temperature case. For example, consider the reduced 
X 2 values for the 2007 fixed surface temperature case in the vicinity of the 
the best-fit values for surface pressure and CH4 mixing ratio (Table HJ). The 
highest reduced x 2 m this subset of parameter space is 2.386, which is lower 
than the 2007 free surface temperature minimum reduced x 2 value of 4.343. 
Thus, at the very least we may claim that the free surface temperature results 
are worse than the one-sigma error bars on the fixed surface temperature 
results. In this case, the disagreement is many sigma. 

In addition to the PGCM parameter sweep of global mean surface pres- 
sure and CH4 mixing ratio, we investigated three different surface radii 
spanning the bounds of pub lished values: 1 180 km (Zalucha et al., 2011a), 



1152 km flElliot et al 



20071 ). and 1130 km (IReinsch et al 



1994J). We found 



that in all cases the 1152 km light curves had the lowest minimum reduced 
X 2 , which differs from the 2D PGCM where Z12 assumed a surface radius of 
1180 km. Table [5] shows the x 2 values for these three radii. 

The best-fit value for CH4 mixing ratio is 1% in all cases fo r fixed surface 
temp erature, which agrees with ZG12 and 



than 



Lellouch et al. 



Young et al 



(119971 ) but is larger 



( 120091 ). The simulated surface pressure in creases from 



1988 to 2006 then drops slightly in 2007. The trend agrees with 



Elliot et al. 
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(12003 



20071), Zalucha et al 



tatively with 



Lellouch et al. 



(2011 a), and ZG12. The values agree quanti- 
d2009|), Zalucha et al. (2011a), and ZG12. We 
point out that our best-fit values for CH 4 mixing ratio are at the upper 
boundary of CH 4 mixing ratio values considered in this study. However, 
the sensitivity of the radiative-conductiv e equilibrium temperature profile 
predicted by the lYelle and Lunind (119891 ) model decreases with increasing 



CH4 mixing ratio. Above about 1%, the model saturates, and little change is 
seen in the magnitude of the radiative-conductive equilibrium temperature 
above this value. Thus, there is little to be gained by running the PGCM at 
CH4 mixing ratio values higher than 1%. 

Figures dHS] show the best-fit modeled and observed light curves. The 
3D PGCM light curves do not contain the oscillations that the 2D PGCM 
exhibited in ZG12, which were due to numerical artifacts in the 2D PGCM 
solution. The models and data agree above the 0.4 normalized flux level. 
At lower flux levels (i.e., closer to the center of the light curve), the model 
first underestimates the flux levels, then overestimates the flux around the 
midtime of the light curve. The difference is approximately 10% or less for 
2002, 2006, and 2007, when noise in the light curve is ignored. The proxim- 
ity of the mismatch to the midtime of the light curve indicates that Pluto's 
lower atmosphere is not being simulated accurately (the model light curves 
probe to a minimum altitude of about 20 km). We may rule out th e uncer- 



tainty in Pluto's surface radius as a cause of discrepancy, because 



zaluchd 



showed that the effect of changing light curve radius was to widen or narrow 
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the model light curve; in fact the flux near the midtime of the light curve 
was affected very little by changing surface radius. The assumption of fixed 
vs. free surface temperature affects the bottom of the light curve only indi- 
rectly through the assumption of energy balance within the column. If the 
discrepancy were due to only haze extinction, then the model would over- 
estimate the flux. Since the model both overestimates and underestimates 
the flux near midtime, the model must instead be somehow misrepresenting 
the temperature gradient. The model flux could be redu ced near the mid- 



time of the light curve if Pluto has a deep troposphere ( jStansberry et al 



19941 ). although it would contradict Zalucha et al. (2011b) who found no 



troposphere. Note that the disagreement in 1988 is larger than the other 



years. The cause of the sharp drop off (kink) at t he 40% flux ley e 



well model ed by both a strong ther mal inversion ( jLellouch et al. 



has been 



extinction ( lElliot and Young 



20091 ) and 



19921 ). but to date no study has conclusively 



proven which scenario is taking place on Pluto. 

3.2. Best- fit PGCM wind, temperature, and surface pressure 

As in the 2D PGCM, the 3D PGCM solution contains essentially no 
vertical motions (the maximum magnitude of vertical velocity in pressure 
coordinates is of order 10" 7 /ibar s _1 ). Similarly, the surface pressure does 
not vary over the globe by more than 10~ 3 /ibar. Unlike the 2D PGCM, 
the meridional motions are non- negligible, occurring in the form of transient 
waves. The maximum magnitude of the instantaneous meridional velocity is 
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2 m s -1 . 

The zonal wind is nearly homogeneous in longitude. Figure [6] shows the 
zonally averaged zonal wind corresponding to the dates of the occultation 
datasets. Three vertical regimes exist in the modeled solution. First is the 
near-surface (below 50 km altitude or above 10 /ibar pressure) with weak 
winds except for a ~2 m s _1 jet centered at 60° latitude. The winds in 
this region are weak because they lie at altitudes subject to surface frictional 
drag. Above this region (between 50 and 130 km altitude and 2 and 10 //bar), 
the winds are easterly (maximum magnitude 6 m s _1 ) between the South 
Pole and 30° latitude and westerly (maximum magnitude 2 m s _1 ) at other 
latitudes. Finally, in the region between 130 km and the base of the Rayleigh 
drag layer at 300 km (0.6 and 2 //bar), the winds are characterized by an 
easterly equatorial jet (maximum magnitude 12 m s -1 ). In the future we 
plan to perform a sensitivity study to determine if the depth of the model 
domain and/or vertical grid resolution affects the vertical boundaries of these 
regions and the associated wind magnitudes. 



2008|) derived an 



Few wind measurements exist for Pluto. ( iPerson et al. 
upper limit of 3 m s _1 , although as explained in ZG12, this analysis was not 
carried out correctly. Moreover, it is unclear what altitude this value applies 
to, making direct comparison difficult. 

The circulation pattern in the 3D PGCM is quite different than in the 2D 
PGCM, which exhibited a zonal circulation in cyclostrophic balance. How- 
ever, we would not expect the 2D and 3D circulations to be the same. The 
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3D PGCM has a day-night insolation contrast which can drive significant cir- 
culations (in addition to the meridional contrast) and also allows for waves 
that can transport energy and momentum (and thus alter the circulation 
pattern) 

Figure [7] shows the zonally averaged temperature corresponding to the 
dates of the occultation datasets. The global temperature field has few hor- 
izontal variations, despite the insolation varying in latitude and longitude. 
The lack of horizontal variability in the temperature field suggests that the 
circulation is efficient at moderating the meridional and day side-night side 
heating gradients. Pluto's rotation period is 6.4 Earth days, while the radia- 
tive relaxation time scale is of order 10 to 100 Earth days. Thus the response 
of the atmosphere is much slower than the daily changes in radiative forc- 
ing, so we would n ot expect a strong day-night temperature contrast to exist. 



Yelle et al. 



( 119951 ) . using assumptions of a frost covered surface in vapor pres- 
sure equilibrium, predicted that the surface temperature and pressure would 
have very little variation on Triton (and Pluto by analogy). This prediction 
is consistent with our GCM results, since the GCM is implicitly adding latent 
heat when the temperature drops below the N 2 freezing temperature and the 
temperature is reset back to the freezing temperature. Future work will focus 
on the contribution of latent heating compared with other heat sources by 
modeling the volatile cycle explicitly. 
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3. 3. Heat transport in Pluto 's atmosphere 



The meridional transport of energy Fe is given by (IPeixoto and Port . 



19921 ) 



[Fe] = c p \Tv\ +g[zv] + l [(u 2 + v 2 )v\ +L[vq], (7) 

where brackets denote a zonal average, overbars denote a time average, v is 
the meridional velocity, g is the gravitational acceleration at the surface, u 
is the zonal velocity, L is the specific latent heat of sublimation, and q is the 
specific humidity. The terms on the right hand side of the equation, from left 
to right, are the transport of sensible heat, potential energy, kinetic energy, 
and latent heat. Since we do not explicitly treat condensables, the latent 
heat term is zero. 

Figure El shows the PGCM values for vertically averaged meridional heat 
transport corresponding to the terms in Eq. [7J At the time of every occulta- 
tion event, the kinetic energy is much smaller than the other terms, which is 
typical of other planetary atmospheres. Pluto's equinox occurred in the year 
1989, and its solstice will occur in 2030. The equations for solar zenith angle 



Peixoto and Port 



19921 ) show 



as a function of latitude and season (see e.g., 
that bodies orbiting the Sun exhibit an equinoctial pattern for a compara- 
tively short fraction of their orbit; the remaining time is spent in a solstitial 
forcing pattern. Thus, by the year 2002, Pluto is already experiencing a sol- 
stitial forcing. For the years 2002, 2006, and 2007, the energy transport is 
always positive. In these years, the South Pole is the summer pole and the 
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amount of diabatic (i.e., solar) heating decreases from south to north. Thus, 
a positive energy transport corresponds to warm air moving northward, or 
a redistribution of energy from south to north. This is consistent with the 
meridionally homogeneous temperature gradient seen in the PGCM temper- 
atures (Fig. [7b-d). Also, in the years 2002-2007, the sensible heat transport 
is larger than the potential energy transport. 

In the year 1988 (Fig. [8^), the total energy is less than in the years 2002- 
2007. The 1988 diabatic heating is equinoctial, with the maximum heating 
near the equator. Thus, we would expect energy transport away from the 
equator, i.e., positive in the Northern Hemisphere and negative in the South- 
ern Hemisphere. Indeed, this behavior occurs in the potential energy term. 
However, the sensible heat transport opposes this motion and is large enough 
to make the total energy positive everywhere. Since the equinox for other 
planetary atmospheres such as Earth or Mars is a transition period, it may 
be that many non-zero transport processes are occurring to nearly cancel 
each other out. A second explanation is that in the PGCM, there is an 
additional energy source when the temperature falls below the freezing tem- 
perature. The temperature is instantaneously snapped back to the freezing 
temperature, implying the addition of latent heat in the atmosphere. We 
suspect that in the 1988 model results, the latent heat term is so strong that 
it is dominating the other terms, and speculate that when an explicit volatile 
cycle is included, it will be important near equinox. 

For v and some variable A (i.e., T, z and {u 2 + v 2 ) in Eq. |7J), we may use 
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the relationship 



= [v] p] + v*A + pA 7 ] , (8) 



where stars denote a departure from the zonal mean and primes denote a de- 
parture from the temporal mean. From left to right, the terms on the right 
hand side of the equation are the transport associated with the mean merid- 
ional circulation, stationary eddies, and transient eddies. Figure [9] shows the 
breakdown of the total energy into the terms of Eq. [HJ In all four years, the 
the contribution from stationary eddies is negligible, which is not unexpected 
since we have included no longitudinal asymmetries in the topography, sur- 
face albedo, emissivity, or surface thermal inertia fields. In 2002 and 2006 
the mean meridional circulation and transient eddies are of comparable mag- 
nitude, while in 2007 the transient term is larger. Again in 1988 there is no 
clear pattern for the energy transport terms. 

3.4. Triton GCM results 

Figure [10] shows the zonally averaged temperature and zonal wind TrGCM 
results for 4 November 1997. Triton's globally averaged temperature is much 
colder (over 30 K) than Pluto and the temperature decreases with height 
rather than increases, due to the assumption that the CH 4 mixing ratio is 
lower by a factor of 50. No te that other radiative-conductive models of Tri- 



ton (e.g. 



Elliot et al. 



20001 ) exhibit temperature increasing with height, as 



does ours when a higher CH4 mixing ratio is used. In the TrG CM assump- 



tions, we have incorporated more recent occultation data from lElliot et al. 
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(119971 ) and spectral data from 



Lellouch et al. 



(120101 ). which results in a lower 



CH4 mixing; ratio . On the other hand, inversion temperatures derived by 



Elliot et al 



(120031 ) of the 1997 HST light curve do not show an inversion in 
the region from 17 to 57 km altitude. Instead, the temperature profiles show 
a complicated vertical structure, but lie in a vary narrow temperature range 
(51-53 K). However, this temperature is still warmer than our results. 

Like in the PGCM result s, the meridional temperature gradient in our 



TrGCM results is weak fas in 



Yelle et al. 



19951 ) . The zonally averaged zonal 



wind contains an easterly jet in the tropics and midlatitudes that extends 
from the top of the assumed troposphere (8 km) to the bottom of the Rayleigh 
friction layer (near 86 km altitude). The jet has a maximum magnitude 
of 4 m s _1 at the equator and around 10 km altitude (10 //bar) pressure. 
The zonal winds at the poles are weakly westerly and achieve a maximum 
magnitude at 70-80 km. Again, the winds in the frictional boundary layer, 
which happens to also extend to 8 km, are weak. 



1989J), de- 



Voyager 2 detected two well-observed plumes (iSmith et al. 
noted as the "west" and "east" plumes, located at —50° and —57° latitude, 
respectively. The west plume appears as a dark column that abruptly ter- 
minates at 8 km altitude. Connected to this plume is a more diffuse clouds 
that extends westward for at least 150 km. Similarly, the east plume also 
rises to an altitude of 8 km and is directed westward at altitude. These 
plumes act as a tracer a nd indicate th a t the flow is easterly, which agrees 



with the TrGCM results. 



Hansen et al. 



( 119901 ) identified a terminator cloud 
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in Voyager 2 images that moved 13m s _1 eastward and was located at an 
altitude of 5 km. The direction and speed of this cloud is opposite in di- 
rection a nd much higher in speed than indicated by the TrGCM results. 



However, 



Hansen et al 



(Il99d ) point out that because the inferred velocity of 
this cloud is so similar to the velocity of the terminator of Triton's surface, 
the cloud may in fact be a stationary, elongated east-west cloud that is being 
illuminated at different points along the cloud by the low altitude sunlight. 
This interpretation would be consistent with the TrGCM results, since in 
the Southern Hemisphere (the observed location of the cloud), the winds are 
ne arly stationary . 



Hansen et al. 



( 119901 ) also observed several crescent streaks at 1-3 km al- 
titude and surface (< 1 km altitude) streaks that potentially indicate wind 
direction at these altitudes. The direction of the crescent streaks was nearly 
uniformly westward, while the surface streaks were variable but mainly north- 
eastward. In the lowest layer of the TrGCM, there is no favored direction 
at the latitudes and longitudes of the crescent and surface streaks; small- to 
medium-scale structures of convergence/divergence and rotation are present. 
Note that the TrGCM is designed to model the full depth of the atmosphere 
and has a surface drag scheme that represents the bulk effect of the frictional 
boundary layer on the free atmosphere. To explain the crescent and surface 
streaks seen by Voyager 2, a mesoscale model the better represents small- 
scale turbulence near the surface is required. Thus, it is not appropriate to 
compare the crescent and surface streaks with the TrGCM results. 
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Like in the PGCM results, in the TrGCM results there is practically no 
vertical motion and the surface pressure v ariation is small c ompared to the 



Yelle et al 



19951 ). Note that the 



globally averaged value (in agreement with 
vertical rising motion observed within Triton's surface plumes is not incon- 
sistent with our result for a number of reasons. First, if the plume material is 
composed of a different material and/or is warmer than the background atmo- 
sphere, this would allow the plumes to be buoyant. Second, the plumes may 



move vertically if they have a non-zero upward initial velocity (see 



Kirk et al. 



1995 



for a detailed discussion of plume sources). Thirdly, our GCM reso- 
lution is much coarser than the horizontal area of the plumes. The GCM 
may be averaging out any upward motion from the plumes with downward 
motion elsewhere in the grid box, resulting in zero net velocity. 

Figure [11] shows the calculated TrGCM light curve calculated and the 
light curve data from the HST on 4 November 1997. The model systemati- 
cally underestimates the data near the midtime by 5-10%, unlike the Pluto 
model light curves where the opposite was the case. Again, the lower atmo- 
sphere is apparently not being simulated accurately. An underestimate by 
the model cannot be accounted for by haze, but could be the result of the 
vertical temperature gradient dropping off too strongly with altitude (i.e., 
the true lower atmospheric vertical temperature gradient is more isothermal 
or contains an inversion). We have investigated the model underestimate by 
performing runs with surface pressure doubled, CH4 mixing ratio increased 
by an order of magnitude, tropopause height doubled, and surface tempera- 
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ture fixed at 52 K (following lElliot et al.l . 120031 ) . None of these modifications 
were successful at increasing the model light curve flux near midtime. 



4. Discussion 

4-1. The effect of Pluto's lapse rate on circulation 

In the 2D and 3D PGCMs (also the TrGCM), we have found practically 
no vertical motion. This is unique versus other planetary atmospheres in 
the Solar System, such as Earth, Venus, Mars, Titan, and the giant plan- 
ets. These planets have one or more overturning circulation patterns such as 
Hadley, Walker, or eddy driven circulations. It is a well-known fact that the 
atmosphere resists vertical motions when a statically stable (subadiabatic) 
vertical temperature gradient is present. In this section, we show that the 
temperature inversion in Pluto's stratosphere appears to be what is prevent- 
ing these vertical motions. 

To demonstrate the effect of a less statically stable atmosphere on Pluto's 



atmos pheric circulation, we replace the non-LTE, shortwave dominated 



Yelle and Lunine 



( 119891 ) radiative-conductive forcing with a simple radiative-convective scheme 
that is representative of a greenhouse dominated atmosphere such as Mars 
or Earth. The external heating term is given by 



^ = -k T (T - T eq ) , (9) 



where kj< is the radiative heating rate and T eq is the equilibrium temperature. 
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T eg is defined by (Zalucha et al., 2010) 

Q [0.5 + 0.75r(p)] r < r t 

-t;, = { (io) 

Q o [0.5 + 0.75T t }[T(p)/T t } 4R/c " r>n, 

where Q Q is the insolation, p is pressure, r is the long wave optical depth, r t 
is the value of r at the top of the convective layer (found by energy balance 
constraints), and R is the specific gas constant. Equation [TU] states that the 
equilibrium temperature follows an adiabat in the surface convective layer 
of the atmosphere and a radiative equilibrium profile above. The long wave 
optical depth is linearly dependent on pressure as r = r (p/p ), where r Q and 
p are reference optical depth and pressure, re spectively. 



In this example (following IZaluchal . 



20101 ) . we assume that the surface 



pressure (set equal to p Q ) is equal to 13.2 /ibar, which is an intermediate value 
for the period 1988-2007. r D = 2 and k? are tuned until Pluto-like temper- 
atures are achieved. The exact values for p a , r , and k^ are not important; 
our main goal is simply to drastically alter the temperature stratification. 

Figure [12] shows a zonally and time averaged latitude-height tempera- 
ture cross-section for Pluto with a greenhouse dominated atmosphere near 
equinox and solstice. The heating profile produces a temperature structure 
that decreases monotonically with height, although some regions with tem- 
perature inversions exist due to dynamical processes. Figure [13] shows the 
corresponding zonally and time averaged latitude-height cross section of mass 
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stream function. There is clearly an overturning circulation pattern (Hadley 
cells); at equinox there are two cells (one cell in each hemisphere) and at 
solstice there is a single planet-wide cell. It is not surprising that with a 
Mars-like temperature structure we obtain a Mars-like circulation. However, 
the preceding example shows the importance of the vertical lapse rate on 
atmospheric circulation. In the nominal PGCM with a steep temperature 
inversion, the mass stream function is effectively zero everywhere. It is ab- 
solutely critical to include the steep temperature inversion in atmospheric 
models of Pluto to ensure a circulation with no vertical motions. 

In models that include explicit surface-atmosphere mass exchange, a lo- 
cal deficit of mass (i.e., low pressure) is created where surface ice deposition 
occurs, and a local abundance of mass (i.e., high pressure) is created when 
sublimation occurs. On Mars, these pressure gradients drive global flows 
from high to low pressure at all altitudes. So-called condensation flows prob- 
ably exist on Pluto, but they will likely be confined to the atmosphere very 
near the surface. The surface exchanges mass with air immediately above 
the surface, and with the radiative-conductively forced vertical temperature 
gradient in place, mass will not be able to move upward or downward to or 
from layers aloft at large scales. Moreover, the weak dayside-nightside tem- 
perature contrast in the PGCM results (due to the slow radiative-conductive 
relaxation rate compared to the rotation rate) suggests the nightside of Pluto 
will have a lower frost deposition rate than if there were no atmospheric heat 
transport. 
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An exception to this scenario is the strength of turbulence and vertical 
mixing in the atmosphere near the surface. If the vertical mixing was able to 
mix sublimating surface volatiles high enough into the atmosphere, then the 
condensation flow need not be confined to the surface. Our GCM does not 
include a parameterization for small-scale turbulence; however its properties 
on Pluto are not likely to known in the near future anyway since it requires 
in situ measurements and empirical modeling to characterize it, which has 
currently only been done for Earth and Mars. 

The PGCM light curves generated from the simulation with fixed surface 
temperature at the freezing point had slightly reduced \ 2 values compared 
with an involatile surface allowed to vary in time and space. A fixed surface 
temperature at the freezing point suggests that Pluto has a large (i.e., thick) 
surface volatile ice reservoir that is able to maintain the freezing temperature 
despite solar and atmospheric forcing. 

4-2. Triton's lower atmosphere 

The TrGCM shows that the lowest ~5 km of Triton's equatorial atmo- 
sphere is characterized by relatively weak zonal winds, before a sharp increase 
in wind speed above this altitude. This behavior is very similar to the Voy- 
ager 2 observation of plum es on Triton that e x tended to 8 km altit ude before 



being blown downstream (ISmith et al 



19891 ) . 



Smith et al. 



( 119891 ) suggested 



that the behavior of the plumes is controlled by a temperature inversion at 
the tropopause and/or vertical structure in the wind speeds. 
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Two important model elements are located in the lowest region of the at- 
mosphere. First, the radiative-conductive-convective heating/cooling scheme 
is set so that in the steady state, the atmosphere below the tropopause (fixed 
at 8 km) follows a moist (with respect to N 2 ) adiabat (i.e., decreases with 
height). Second, a frictional layer is present from the surface to a pressure 
level equal to 70% of the surface pressure (in this case 17.64 //bar). When 
full model dynamics and radiative-conductive-convective heating/cooling are 
allowed, we find that the temperature actually increases slightly with height 
in the lower atmosphere. To test whether it is the external heating/cooling 
or surface friction that is dominating the lower atmosphere, we performed 
a simulation where the diffusion of potential temperature was turned off in 
the radiative-conductive-convective heating/cooling scheme. This action re- 
moves the thermally-induced troposphere. In this case the near-surface winds 
were still weak. Thus, the behavior of the plumes once they have been re- 
leased into the atmosphere can be caused by strong wind shear capping the 
lower atmosphere, rather than thermal effects. This resu lt correspond s to th e 



second intepretation of the plume behavior su ggested by 



Smith et al. 



Ingersolll ( 



(1989) 



1990) used Ekman 



At the time of the Voyager 2 encounter, 
layer theory (near-surface balance between the pressure gradient force, Cori- 
olis force, and turbulent drag) to predict an anticyclone at the south pole. 
Anticyclonic flow corresponds to flow opposite to the body's rotation, or in 
this case westerly. Our TrGCM results do show westerly flow at the south 
pole (poleward of ~60°), but the wind speed is less (< 0.5 m s -1 ) than the 
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Ingersolll ( 119901 ) estimate of 5 to 15 m s l . Dust devils, which are not resolv- 
able in our model, hav e also been proposed as an other dynamical mechanism 



to explain the plumes ( llngersoll and Tryka 



1990). However, the static stabil- 



ity of the lower atmosphere of our model would prohibit dust devils, which 
can only form in the presence of neutral or unstable vertical temperature 
profiles. 



5. Conclusion 



We have simulated the atmospheres of Pluto an d Triton using a 3D GCM 



Yelle and Luninel (119891 ) 



based on the MIT GCM dynamical core and the 
and Zalucha et al. (2011b) radiative-conductive-convective schemes. We 
find that, like the 2D version of this GCM, there is no significant large- 
scale vertical motion and hence no overturning cells in the Pluto or Triton 
atmospheres. Triton's plumes, which are not explicitly present in our model, 
can be an exception to this rule if they are composed of a different and/or 
more buoyant material than the background atmosphere or escape the surface 
with a non-zero vertical velocity. 

Some meridional motion is present in our simulations (maximum mag- 
nitude 2 m s -1 ) in the form of transient waves and mean circulations that 
transport heat meridionally very efficiently As a result, the meridional tem- 
perature gradient is weak (as is the dayside-nightside temperature contrast). 
On both Pluto and Triton, a zonal jet is located at the equator and mid- 
latitudes. Pluto exhibits a more complicated vertical structure in the zonal 
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winds with three distinct regimes present. 

The strong temperature inversion in Pluto's lower atmosphere prevents 
large-scale overturning cells. By substantially changing the lapse rate so that 
it was similar in vertical structure to a greenhouse-dominated atmosphere 
(e.g., Mars), we were able to recover Hadley cells. We also showed that the 
macroscopic behavior of plumes on Triton can be explained by the dynami- 
cal properties (i.e., strong wind shear) of the lower atmosphere alone, since 
a layer of abruptly strong winds overlayed a surface layer of essentially qui- 
escent winds when the convective component of the heating external scheme 
was removed from the TrGCM. 

The PGCM and TrGCM are continually improving. In the future we 
plan to include CH 4 transport, a multi-layer surface model, and a frost cy- 
cle. The latter two items represent particular challenges because the sur- 
face and subsurface parameters (e.g., albedo, surface thermal inertia, and 



emissivity) r epresent a particu l arly 



constrained ( Hansen and Paige 



1996 



arge paramet 



Young 



er space that is not well 



2012). These may also require 



multi-year simulations to equilibrate the surface frost thickness and/or sub- 
surface temperature (heat storage). 

Because Pluto and Triton are so similar in their bulk properties it is 
useful to study them in tandem. A relatively large amount of data has 
been obtained through observations of stellar occultations by Pluto, and 
more are planned to be acquired by the New Horizons spacecraft, currently 
en route to Pluto. While Triton has already been visited by Voyager 2, 
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those observations are now over two decades old and Triton's atmosphere has 
potentially undergone seasonal changes. With further observations of Triton, 
such as stellar occultations, spectra, and surface albedo, we will be able to 
better constrain the TrGCM, and have the secondary benefit of improving 
the PGCM. 
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Table 1: Pluto and Triton GCM parameters 



Parameter 


Pluto value 


Triton value 


Surface radius (km) 


1152 a 


1353 


Orbital eccentricity 


0.251 


0.0097 6 


Semimajor axis (AU) 
Rotation rate (1(T 5 s _1 ) 


39.8 


29.9 6 


-1.13856 


-1.2374 


Obliquity with respect to Sun (°) 


60.4 


28.3 


Surface gravitational acceleration (m s _1 ) 


0.63 


0.77 


Wavelength-integrated solar constant (W m~ 2 ) 


0.864 


1.137 


Ecliptic longitude of perihelion (°) 


186 


38.0 



different from ZG12. 



b of Neptune with respect to Sun. 



Table 2: Best-fit PGCM parameters (surface temperature free to evolve) 



Event 


Global mean surface 


CH 4 mixing 


Minimum 


Degrees of 




pressure (//bar) 


ratio (IO -3 ) 


reduced \ 2 


freedom 


1988 KAO 


12±4 


6±4 


6.863 


998 


2002 UH 2.2m 


24±4 


> io 


2.276 


2390 


2006 AAT 


24±4 


> io 


6.773 


3598 


2006 Siding Spring 


24±24 


> io 


1.438 


1798 


2007 Mt. John 


> 28 


6±4 


4.343 


1999 
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Table 3: Best-fit PGCM parameters (surface temperature fixed) 



Event 


Global mean surface 


CH4 mixing 


Minimum 


Degrees of 




pressure (//bar) 


ratio (10 -3 ) 


reduced x 2 


freedom 


1988 KAO 


8 ± 4 


> 10 


5.336 


998 


2002 UH 2.2m 


20± 4 


> 10 


1.994 


2390 


2006 AAT 


20± 4 


> 10 


5.259 


3598 


2006 Siding Spring 


20± 4 


> 10 


1.338 


1798 


2007 Mt. John 


16± 4 


> 10 


2.210 


1999 



Tabic 4: Selected reduced \ 2 f° r the 2007 fixed surface temperature case 



Global mean surface 


Reduced x 2 f° r 


Reduced x 2 


pressure (/xbar) 


CH 4 mixing 


CH 4 mixing 




ratio of 0.6% 


ratio of 1% 


12 


2.386 


2.217 


16 


2.332 


2.210 


20 


2.332 


2.255 



Tabic 5: Minimum reduced x 2 f° r selected surface radii (surface temperature 
fixed) 



Event 


1130 km 


1152 km 


1180 km 


Degrees of freedom 


1988 KAO 


6.492 


5.336 


5.367 


998 


2002 UH 2.2m 


2.000 


1.994 


2.379 


2390 


2006 AAT 


6.046 


5.259 


7.142 


3598 


2006 Siding Spring 


1.391 


1.338 


1.648 


1798 


2007 Mt. John 


2.269 


2.210 


2.384 


1999 
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Figure 1: Top panel: Best PGCM fit for surface temperature free (red curve) and surface 
temperature fixed (blue curve) to f988 KAO normalized light-curve data (black points). 
The dashed horizontal line is the zero flux level. Bottom panel: residual between models 
and data. 
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Figure 2: Top panel: Best PGCM fit for surface temperature free (red curve) and surface 
temperature fixed (blue curve) to 2002 UH 2.2m normalized light-curve data (black points). 
The dashed horizontal line is the zero flux level. Bottom panel: residual between data and 
models. 
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Figure 3: Top panel: Best PGCM fit for surface temperature free (red curve) and surface 
temperature fixed (blue curve) to 2006 Siding Spring normalized light-curve data (black 
points). The horizontal vertical line is the zero flux level. Bottom panel: residual between 
data and models. 
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Figure 4: Top panel: Best PGCM fit for surface temperature free (red curve) and surface 
temperature fixed (blue curve) to 2006 A AT normalized light-curve data (black points). 
The dashed horizontal line is the zero flux level. Bottom panel: residual between data and 
models. 
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1.2 
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Figure 5: Top panel: Best PGCM fit for surface temperature free (red curve) and sur- 
face temperature fixed (blue curve) to 2007 Mt. John normalized light-curve data (black 
points). The dashed horizontal line is the zero flux level. Bottom panel: residual between 
data and models. 



50 




-90 -60 -30 30 60 -90 -60 -30 30 60 



Latitude (degrees) Latitude (degrees) 

Figure 6: Zonally averaged zonal wind (m s _1 ) corresponding to best-fit PGCM results 
from (a) 9 June 1988 (b) 21 August 2002 (c) 12 June 2006 (d) 31 July 2007. 
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Figure 7: Zonally averaged temperature (K) corresponding to best-fit PGCM results from 
(a) 9 June 1988 (b) 21 August 2002 (c) 12 June 2006 (d) 31 July 2007. 
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Figure 8: Simulated vertically averaged meridional transport of energy on Pluto corre- 
sponding to (a) 9 June 1988 (b) 21 August 2002 (c) 12 June 2006 (d) 31 July 2007. The 
black solid line is the total energy transport, the red dashed line is the potential energy 
transport, the green dotted line is the sensible heat transport, and the blue dashed-dotted 
line is the kinetic energy transport. These terms are defined in Eq. [7] and have been mul- 
tiplied by a factor of 2-n a cos <pp s (where a is the surface radius, 4> is latitude, and p s is 
surface pressure) to obtain units of Watts. The results have been time averaged over 90 
Earth days. Positive values indicate northward transport of heat. 
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Figure 9: Simulated vertically averaged meridional transport of energy on Pluto corre- 
sponding to (a) 9 June 1988 (b) 21 August 2002 (c) 12 June 2006 (d) 31 July 2007. The 
black solid line is the total energy transport, the red dashed line is the mean meridional 
transport, the green dotted line is the stationary eddy transport, and the blue dashed- 
dotted line is the transient eddy transport. These terms are defined in Eq. [8] and have 
been multiplied by a factor of 2ira cos <pp s (where a is the surface radius, 4> is latitude, and 
p s is surface pressure) to obtain units of Watts. The results have been time averaged over 
90 Earth days. Positive values indicate northward transport of heat. 
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Figure 10: TrGCM results from 4 November 1997. (a) Zonally averaged temperature (K) 
and (b) zonally averaged zonal wind (m s" 1 ). 



55 



1.2 




200 300 400 500 600 700 800 
Time (s) after 05:00:00 UTon 4 November 1997 



Figure 11: Top panel: light curve corresponding to the TrGCM (blue curve) and HST 
data (red points) on 4 November 1997. The dashed horizontal line is the zero flux level. 
Bottom panel: residual between data and model. 
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Figure 12: Zonally and time averaged temperature (K) for the PGCM with greenhouse 
gas absorber (a) equinox-type case (b) solstice-type case. The surface pressure for this 
example is 13.2 fjbar. 
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Figure 13: Zonally and time averaged mass stream function (10 7 kg s _1 ) for the PGCM 
with greenhouse gas absorber (a) equinox-type case (b) solstice-type case. Positive flow is 
counter-clockwise. The surface pressure for this example is 13.2 /ubar. 
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